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ABSTRACT 

O ■ We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treat- 

jf-H ■ ment of star formation, supernovae feedback, and chemical enrichment, to perform 

c/3 \ eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of 

the progenitors, their rotation axes, and their orbital parameters and analyze the kine- 
matic, structural, and chemical properties of the remnants. Six of these simulations 
result in the formation of a merger remnant with a disc morphology as a result of the 
large gas-fraction of the remnants. We show that stars formed during the merger (a 
sudden starburst occur in our simulation and last for 0.2 — 0.3 Gyr) and those formed 
after the merger have different kinematical and chemical properties. The first ones are 
located in thick disc or the halo. They are partially supported by velocity dispersion 
and have high [a/Fe] ratios even at metallicities as high as [Fe/H] = —0.5. The former 
ones - the young component - are located in a thin disc rotationally supported and 
have lower [a/Fe] ratios. The difference in the rotational support of both components 
results in the rotation of the thick disc lagging that of the thin disc by as much as a 
factor of two, as recently observed. We find that, while the kinematic and structural 
properties of the merger remnant depends strongly upon the orbital parameters of 
the mergers, there is a remarkable uniformity in the chemical properties of the merg- 
■ ers. This suggests that general conclusions about the chemical signature of gas-rich 

mergers can be drawn. 
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1 INTRODUCTION 
1.1 Disc-Disc Mergers 

In a A Cold Dark Matter Universe (ACDM), the forma- 
tion of structures generally involves the merging of smaller 
structures. Mass is built up predominantly by the merger of 
objects of masses between 3 — 30% of t he galaxy's total mass 
at redshift z — ([Stewart et al.ll2008h . Since z = 2, approx- 
imately 60% of galaxies have experienced a major merger, 
which we will take in this paper as those involving a ra- 
tio of masses 3:1 or higher. Also, an important fraction of 
the baryonic mass (30 - 50%) of z — galaxies is provided 
by major mergers (jStewart et al.l l2009). Recent observations 
show that nearly 20% of the massive disc galaxies have un- 



dergone a major merger since z ~ 1.5 (|L6pez-Saniuan et alJ 
2009) This process is virtually scale- free, meaning that low- 
mass dark matter halos, in which disc galaxies presumably 
reside, have undergone similar merging histories to high- 
mass dark matter halos in which, predominantly, ellipticals 
reside. This at first is puzzling, as mergers have long been 
associate d with forming; spheroidal structures, i.e. early-type 
galaxies ([Toomrell 19771 ). Following this line of reasoning, ev- 
ery galaxy in the nearby universe should be an elliptical 
or an irregular, which is clearly not the case. Furthermore, 
observed elliptical galaxies can be described by simple rela- 
tions relating colors, luminosities, masses, etc. These simple 
relations are difficult to reproduce wit h a scenario involving 
collisions (jBekki fe Shiovall 19981 . [l999h . 

The resolution of this, of course, lies with the baryonic 
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physic s. Studies using two stable discs as initial condi- 
tions (ISpringel fe Hernquistl 120051 : iRobertson et al.1 I200& 
Brook et all boO/t hereafter B07: IRobertson fc Bullockl 
2008 ) as well as fully cosmological simulations 
(Go vernato et all 120081 ) have shown that binary gas- 
rich mergers can result in disc morphologies. Hence, major 
mergers between gas-rich disc galaxies can possibly play 
an important role in the formation of a large number of 
the spiral galaxies we observe today. These facts bring an 
addition to the original scenario for the formation of disc 
galaxies, which involved the assembly of a halo of dark mat- 
ter and ionized gas by hierarchical clustering and mergers of 
sub-halos, followed by the collapse and virialization of the 
dark matter halo, and finally the dissipative collapse of the 
baryonic component, resulting in a rotation ally-supported 
disc, that will later fragment to form stars ([White Sz Reesl 
1 19781 : iFall fe Efstathioulfl980l : iBlumenthal et al]ll984h ~ 

One must also consider a formation scenario involv- 
ing major collisions between spiral galaxies. The first nu- 
merical simulations of major collisions between disc galax- 
ies were performed under the assumption that dissipa- 
tion by gasdynamical effects was negligible. The results 
showed that the remaining galaxy had several properties at- 
tributed to elliptical ga l axies, including luminosity profiles 
( Naa b fe Burkertl 120031 : iBarned l200i IBarnes fc Hernquistl 
1992 : iHernauistl Il99l Il993l: iBekki & Shioval [l998l . Il999l : 
Naab, Burkert, &; Hernquistl 1 19991 ). Even with the quite 
recent inclusion of feedback processes, the results re- 



mained mostly unchanged (IBarnes &; Hernquist J 1996; 
iMihos fc Hernquistlll996l : ISpringelll2000l : iNaab et alJl2006h . 
Minor collisions have also been studied in the past 
and results showed that this kind of collision tends to 
form galaxies whose disc is destroyed or very perturbed 
(jHernquist fe Mihosl Il995l : iBekkil Il998l ). All the aforemen- 
tioned studies assumed a ratio gas/stars more representa- 
tive of present low redshift mergers. Of course, stars have 
not always been present inside galaxies. Going back in time, 
we can expect to find fewer and fewe r stars, and therefor e a 
more important gaseous component ([Stewart et al.l 2009). 

1.2 Thick Discs 

The Milky Way's thick disc is made predominantly of 
old stars, with ages in the range of 8 — 12Gyr, signifi- 
cantly older than the typical star found in the thin disc 
(|Gilmore. Wvse. fe Jonesl Il995h . The thick and thin disc 
differ in their structural, chemical, and kinematical prop- 
erties: (1) The thick disc has a scale heigh t of about 
900 p c, compared to 300 pc for the thin disc ([juric et al.l 
2008). (2) The stars in the thick disc have a low metallic- 
ity, [Fe/H] ~ —0.6 on average, and a ratio [a/Fe] larger 
than the stars in the thin disc. (3) Stars in the thick disc 
tend to ro tate slower than stars in the thi n disc, by about 
40 km s _1 ([Gilmore, Wvse, &; Kuiikenll 19891 ) and have an av- 
erage vertical velocity dispersion of 45kms~ 1 compared to 
~ 10 — 20kms — 1 for the stars in the thin disc (Delhave 1965; 
IChiba fe Beersll2000h . 

The formation of the Milky Way's thick disc re- 
mains an open question, with popular theories in- 
cludi ng the heating of a thin disc by minor merg- 
ers ( Quinn, Hernquist. fe Fullagarl Il993l: Kazantzidis et al. 



2008: ]Kazantzidis. Zentner. & Bullock 2008), or the direct 



accretion of stars, which are tidally torqued into the plane 
of the existing thin disc (|Abadi et al.l l2003h . or d ynamical 
heating by destruction of very massiv e star clusters ([Kroupal 
120021 : lElmegreen fc Efmegreenl 12006). We favor a scenario 
in which the Milky Way's thick disc is formed during the 
high-redshift epoch in which mergers are most common in 
ACDM, and accretion rates are high, with the disc which 
formed during this turbulent period being born thick; the 
thin disc subsequently grows during the relatively quiescent 
period which follows, which in th e case of the Milky Way is 
appr oximately the last 10 Gyrs (|Brook et al.1 l2004bL 120051 . 
120061 ). 

Recent observations suggest that all spiral galaxies are 
surrounded by a red, flattened envelope of stars, and these 
structures have been interpreted as s howing that thick discs 
are u biquitous in spiral galaxies ([Dalcanton &; Bernstein! 
2002). Hence, understanding the formation of the thick disc 
will provide major insights into the formation processes of 
spiral galaxies. If we make the hypothesis that thick disc are 
formed during gas-rich major mergers, it would mean that 
a significant number of spiral galaxies have experienced ma- 
jor mergers. Indeed, simulations have sho wn that gas-rich 
mergers tend to form two disc components ([Robertson et al.l 
2006; B07). These simulations included the effect of star for- 
mation and supernovae feedback, which allowed to study the 
kinematical and structural properties of galaxy merger rem- 
nants in detail. However, they did not include a detailed 
treatment of chemical enrichment and, therefore, could not 
predict the abundances of metals in stars and the inte rstel- 
lar medium (but see also IBekki Shioval [T998. 1999, who 
studied the formation of ellipticals galaxies due to a more 
violent major merger) . It is well-established that the thick 
disc of the Milky Way possesses a unique chemical signa- 
ture, which is believed to be directly related to the forma- 
tion process and it is important to analyze if the merger 
scenario for the formation of the thick disc component is 
compatible with these chemical signatures. This is the basis 
of our study here: we believe that the remnants of major 
mergers may be associated with a significant fraction of th e 
"red envelopes" observed in iDalcanton &; Bernstein! (|2002h , 
i.e. with extragalactic thick discs. Since we expect the con- 
ditions for merging galaxies during hierarchical clustering at 
a high redshift to vary, we need to consider a wide range of 
parameters, such as mass ratio, orbit, and gas fraction of 
the merger progenitors. In this paper, we present a series 
of 8 simulations of mergers of gas-rich disc galaxies, with 
different mass ratios and orbital parameters. We analyze 
the kinematic and structural properties of the merger rem- 
nants. The results for one of the simulations have already 
been presented in a previous paper (B07). We now complete 
this work by presenting the entire suite of simulations. 

For kinematic properties, we determine the relative im- 
portance of rotation and velocity dispersion in providing 
support. We also investigate if a counter-rotating component 
arises naturally in a ga s-rich merger, as such co mponents are 
observed in nature (|Yoachim &; Dalcantodl2006h . For struc- 
tural properties, we determine if the luminosity profiles of 
these thick discs formed in major mergers follow an expo- 
nential law, and if so, what their respective scale lengths are. 
We determine if the thick and thin discs in merger remnants 
are coplanar, or if there is a significant angle between them, 
something that was not included in our previous work. 
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However, the main goal of the present paper is to study 
the viability of the merger mechanism for the formation of 
the thick disc from the chemical evolution's point of view, 
but without neglecting the kinematical and dynamical prop- 
erties. With this in mind, we choose initial conditions likely 
to produce a discy remnant. 

In this paper, we use the numerical algorithm GCD+, 
which includes a detailed treatment of chemical enrichment, 
to investigate whether key abundance ratios such as [a/Fe] 
and [Fe/H] in simulated remnants of gas-rich disc galaxy 
mergers have values and distributions similar to the ones 
found in the observed disc galaxies, including the Milky Way. 

The reminder of this paper is organized as follows: in 
§2, we briefly describe the numerical simulations employed, 
including the N-body/SPH code adopted (GCD+), the ini- 
tial conditions software (GalactlCS), and the basic param- 
eters of our 8 realizations. The kinematical and structural 
properties of the merger remnants are presented in §3 and 
discussed in §4. Conclusions are presented in §5. 



2 THE NUMERICAL SIMULATIONS 
2.1 The Algorithm 

All simulations were performed using GCD+ 
(jKawata fe Gibson! l20Q3al lbh. self-consistently modeling 
the effects of gravity, gas dynamics, radiative cooling, 
and star formation. GCD+ is a Tree/SPH algorithm that 
includes Type la and Type II SNe feedback, and traces 
the lifetimes of individual stars, enabling us to monitor 
the chemical enrichment history of our simulated galax- 
ies. Star formation occurs in a convergent gas velocity 
field where gas density is greater than a critical density, 
p crit = 2 x 10 _25 gcm -3 . The star formation rate (SFR) 
of eligible gas particles is then dp*/dt = dp g /dt = c*p g /t g 
where c* = 0.05 is a dimensionless star formation efficiency, 
and t g is the dynamical time. This formula corresponds 
to a Schmidt law: SFR oc p 1 - 5 . The mass, energy, and 
heavy elements are smoothed over the neighboring gas 
particles using the SPH smoothing kernel. The code 
separately tracks the abundances of 9 elements: H, He, C, 
N, O, Ne, Mg, Si, and Fe. Gas within the SPH smooth- 
ing kernel of Type II SNe explosions is prevented from 
cooling, creating an adiabatic phase for gas heated by 
such SNe. This adiabatic phase is assumed to last for the 
lifetime of the lowest mass star that ends as a Type II 
SNe, i.e., the lifetime of an 8M star (100 Myr). This 
i s sim ilar to a model presented in iThacker fe Couchmanl 
(2000). Stellar yields are dependent of the progenitor 
mass and metallicity. For Type II SNe, we use the 
metallicity- dependent stellar yields of IWooslev fe Weaver] 
(1995) for stars with mass over 11M©, while for low- 
and intermediate-mass stars, we use the stellar yields of 
Ivan den Hoek fe Groenewegenl (Il997h. For Type la SNe, 



ad opt the model of Kobavashi. Tsuiimoto. fe Nomotol 
(2000) in which Type la SNe occur in binary systems with a 
primary star and a companion of a definite mass and metal- 
licity range. The yields are from llwamoto et al.l (|l999l ). Note 
that GCD+ relaxes the instantaneous recycling approxima- 
tion and takes into account the lifetime of progenitor stars, 
and chemical enrichment from intermediate mass stars. The 



code calculates, within each timestep, the amount of energy 
and heavy element released by every star particle and 
distributes them over the neighbours using the SPH kernel 
for both Type II and la SNe. We did not include in the 
present paper any black hole feedback. Several studies have 
demonstrated that this can be very important for the final 
morphology and th e termination of the star formation in the 
merger remnant (Springel, Pi Matteo fe Hernauistl 120051 : 
lOkamoto et al.l [20051 : 1 Johansson et al.l 120091 ) . However, the 
implementation of the blackhole feedback in numerical 
simulations is still highly ambiguous and we prefer here, for 
simplicity, to include only SN feedback. For more details 
about GC D+, we refer the reader to iKawata fe Gibsonl 
(|2003al lbh: lBrook et all (|2004ah . 



2.2 Setting Up Initial Conditions 

In each simulation, the initial conditions consist of two 
galaxies with exponential gas discs embedded in dark mat- 
ter halos. The se galaxies are generated using the Galac- 
tlCS package ([Kuiiken fe Dubinskil Il995l ) , producing sta- 
ble discs for a large number of galactic rotations. For the 
dark matter halos, GalactlCS use the lowered Evans model 
(Evansl [l993l ). which leads to a constant-density core. This 
model differs from models which have a steeper central den- 
sity profile suggested by ACDM cosmological simulations 
dNavarro. Frenk. fe Whrtel Il996l. Il997l : iMoore et alJ Il999l : 
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iGhigna et al.l l2000l : Ijing fe Sutol l200of iKlvpin . et^l200ll ). 
For our purpose this difference is not critical, as the cen- 
tral regions contain only a small fraction of the total mass. 
The star formation mechanism is shut down for the first 
~ 180Myr enabling very gas-rich mergers. In our simula- 
tions, the gas ratio immediately before merger is more rel- 
evant than the initial value. Pure gas disc will induce an 
unrealistic starburst initially. However, by turning off star 
formation initially, we can control the gas fraction at the 
merger epoch. For each simulation, the most massive galaxy 
(Gall) has total mass of 5 x lO n M0 (except for simula- 
tion Mil) and a scale length of 4.5 kpc. The mass and scale 
length of the second galaxy (Gal2) depends on the specific 
mass ratio chosen for the simul ation. Using the observed 
gas fraction of g alaxies at high (lErb et alj_ _2006) and low 
(jMcGaughl 120051 ) redshift. Istewart et all (|2009h provide a 
relation between gas fraction and redshift for galaxies of 
various masses. Using this relation (their equation [1]) al- 
lows us to estimate that the mergers in the epoch between 
z = 2.37 — 3.7 would be between galaxies which have gas 
fractions as high as the ones used in our study. Note that 
this is an upper bound of su ch epoch due to t he fact that the 
high redshift observations of lErb et al.l (|2006h are considered 
lower bounds of the gas content. Also, these redshifts coin- 
cide with the epoch of formation of the stars in the Milky 
Way thick disc. 

The initial configuration of the system is illustrated 
in Fig. [1] Gall is initially located at Ri = (50, 0, 0) kpc, 
with no initial velocity, and its rotation axis pointing in 
the Z direction. Gal2 is located at R2 = (0, 0, 15)kpc (ex- 
cept for simulations M12C and M12orb), with initial velocity 
V2 = (0, 100, 0)kms _1 (except for simulation M12z), and its 
rotation axis in the X — Z plane at an angle relative to the 
Z-axis. We shall refer to this angle as the interaction angle. 
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Figure 1. Geometry of the initial conditions. The y-axis is pointing away from the viewer. Gall is initially at rest, while Gal2 is initially 
moving away from the viewer, in the -\-Y direction, as the thick arrow shows. The dotted lines indicate the rotation axes of the galaxies. 
Gall is in the X — Y plane, with the rotation axis in the Z-direction. The rotation axis of Gal2 is in the X — Z plane, at an angle 
relative to the Z-axis. Both galaxies are rotating clockwise when seen from above, hence the left edges are moving away from the viewer, 
and the right edges are moving toward the viewer. The particular case shown on this figure corresponds to simulation M12, with the 
edges of the discs located at two scale-lengths. 



All encounters are prograde-prograde, except for simulation 
rM12, for which the rotation of Gal2 is retrograde. 



2.3 A Series of 8 Gas-Rich Mergers 

We have performed a total of 8 simulations with different 
initial conditions, in which each galaxy is modeled using 
100,000 dark matter halo particles and 40,000 gas particles, 
for a total of 280,000 particle s. This resolution is low er by a 
factor 2 than the one used in Robertson et al . (2006). bmce 
the primary goal of this work is to experiment on chemical 
enrichment following a major merger, we had to make some 
choices relatives to the numerical resolution we had to use. 
The primary concern being the relatively high computation 
time needed to distribute enrich gas to the particles in the 
neighborhood. The baryon/DM mass fraction is 17%, which 
is equal t o the universal rati o Ob /Op according to recent 
estimates (jSpergel et al.ll2007l ; iKomatsu et al.ll20Q9h . Every 
simulation starts with an initial metallicity log(Z/Z ) = — 4 
and a-element abundance [a/Fe] = 0.30. 

Table [1] lists the parameters used for generating the 
initial conditions for each simulation. The five simulations 
within the M12 "family" involve a major merger of two 
galaxies with a mass ratio of 2:1. We start with M12, which 
has already been discussed in detail in B07. The galax- 
ies collide after 320 Myr, and by that time 6% of the gas 
has already been converted to stars, leaving a gas frac- 
tion /gas = 0.94. The next 4 simulations within that fam- 
ily are variations of the "parent" simulation M12. In sim- 
ulation M12orb, we increase the initial separation between 
the galaxies, and in simulation M12z, we add an additional 
velocity component in the z-direction. In both cases, the 
merger is delayed, and this results in lower gas fractions at 
the time of the merger. In simulation Ml 290 we changed the 
interaction angle to 90 degrees. Simulation rM12 is similar 



to simulation M12, except that the direction of rotation of 
Gal2 is reversed, making it retrograde relative to the orbit. 

With the last 3 simulations, we consider different mass 
ratios. In simulation Mil, we used equal-mass galaxies. Also, 
this simulation is the only one for which we use a smaller 
mass for Gall: 2.5 x 10 n M Q instead of 5 x 10 n M Q . In simu- 
lation M13, we use a mass ratio of 3:1. This is still considered 
to be a major merger. In the simulation Ml 10, we use a mass 
ratio of 10:1, which corresponds to a minor merger. 

Column 8 of Table 1 gives the initial eccentricities of 
the orbits. Except for simulation Ml 10 (the minor merger), 
the eccentricities are quite high, corresponding to elongated 
elliptical orbits. As a result, the pericenters are smaller than 
the actual size of the galaxies, and the collision takes place 
before a full orbit is completed. Column 9 gives the angle be- 
tween the position of the pericenter and the line of nodes (in- 
tersection b etween the orbital plane a nd the plane of Gall; 
see p. 632 of lToomre fe Toomrel[l972h . Column 10 gives the 
inclination i of Gal2 relative to the plane of the orbit. 

For the simulation Mil, the softening length of the al- 
gorithm is 0.65 kpc for the dark matter and 0.48 kpc for the 
baryons. For the other simulations, the softening length is 
0.82 kpc for the dark matter and 0.61 kpc for the baryons. 

Table [2] lists the time £/ elapsed between the begin- 
ning and the end of the simulation, the time t co ii elapsed 
between the beginning of the simulation and the collision, 
and the final gas fraction /g^ al . In each simulation, the col- 
lision leaves a very complex and chaotic system, which then 
relaxes to a "quiescent state" where the final structure and 
kinematics are well-established, and the star formation rate 
has dropped significantly. We stopped the simulations after 
2.5 Gyrs, when the quiescent state had been reached. One 
simulation (M12orb) was extended to tf = 3.40 Gyrs, which 
enabled us to check that the final state was indeed quies- 
cent. Comparing the values of fg^ 1 with the values of / gas 
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Table 1. Initials Conditions for All Simulations 



Run 


M Ga u(M ) 


mass ratio 


/gas 


e 


V 2z (kms- 1 ) 


R 2 (kpc) 


e 


a; [degree] 


/.[degree] 


M12 


5.0 x 10 11 


2:1 


0.94 


30.0 





0, 0, 15 


0.79 


90 


-13 


M12orb 


5.0 x 10 11 


2:1 


0.67 


30.0 





0, 50, 15 


0.86 


51 


-13 


M12z 


5.0 x 10 11 


2:1 


0.83 


30.0 


100 


0, 0, 15 


0.66 


29 


-20 


M1290 


5.0 x 10 11 


2:1 


0.96 


90.0 





0, 0, 15 


0.91 


90 


73 


rM12 


5.0 x 10 11 


2:1 


0.91 


210.0 





0, 0, 15 


0.79 


90 


167 


Mil 


2.5 x 10 11 


1:1 


0.88 


30.0 





0, 0, 15 


0.89 


90 


-13 


M13 


5.0 x 10 11 


3:1 


0.92 


30.0 





0, 0, 15 


0.70 


90 


-13 


MHO 


5.0 x 10 11 


10:1 


0.85 


30.0 





0, 0, 15 


0.37 


90 


-14 



Table 2. Lookback Times and Final Gas Fraction 



Run 


t f [Gyrs] 


*coii [Gyrs] 


j? final 
J gas 


M12 


2.50 


0.27 


0.14 


M12orb 


3.40 


0.83 


0.18 


M12z 


2.50 


0.50 


0.20 


M1290 


2.50 


0.25 


0.10 


rM12 


2.47 


0.32 


0.18 


Mil 


2.50 


0.40 


0.06 


M13 


2.50 


0.32 


0.18 


MHO 


2.50 


0.42 


0.30 



listed in Table \T\ reveals the efficiency of the starburst in 
converting gas to stars. The gas fraction typically dropped 
from ~ 90% to ~ 15%. It is interesting to note that simula- 
tion Ml 10, the only case of a minor merger, was less efficient 
in converting gas to stars than the other simulations, with 
/ g fi a n s al = 0.297. 



3 RESULTS 

Following the approach of iRobertson et all (2006), we have 
computed the SFR for all the simulations, and used these 
results to identify stellar populations. Some examples of 
them are shown in the top left panels of Figs. [BHSJ No- 
tice that the simulations start at t = 0. The results are 
quite similar for all simulations. A major starburst invari- 
ably occurs during the merger. This starburst peaks between 
~ 6OM0 yr _1 and 500M©yr -1 , compar able to that see n 
in high-redshift Lyman-break galaxies (|Erb et al.1 l2QQ6h . 
The results are qualitatively similar to recent work involv- 
i ng similar gas ratio but different feedback assumptions 
([Johansson et l2QQ9h o r different gas ratio and different 
feedback (|Cox et al.ll2008l ) . The main quantitative difference 
come from t he amount of ga s avail able during the merger. 
As stated in I Johansson et all ([20091 see their Figs. 1 and 3), 
there is a huge difference in the maximum value of the SFR 
depending on gas mass fraction. The SFR is also dependent 
of the numeri cal methods us ed for feedback in case of ma- 
jor mergers ([Cox et alJl2006h . After the starburst, the star 
formation steadily drops, as the system relaxes and a disc 
forms. We identify the beginning and the end of the merger 
with the beginning and the end of the starburst, respectively. 
We then define two populations of stars: old stars, which in- 
clude stars already present in the galaxies before the merger 
and stars formed during the merger by the starburst, and 
young stars, which include all stars formed after the merger, 



when the starburst is completed. The definitions for young 
and old stars remain the same throughout the paper. The 
dashed vertical lines in the top- left panels of Figs. [BH9] indi- 
cate the beginning and end of the starburst, identified by eye 
as the merger boundaries are not crucial for the subsequent 
analysis. The time when the starburst begins is listed in Ta- 
ble [2] For most simulations, the starburst lasts 0.2-0.3 Gyr. 
The simulation M12z (Fig. is a notable exception, with a 
starburst that is weaker and lasts for 0.5 Gyr. 

We use the term remnant to designate the final state of 
each simulation. Our simulations are meant to represent the 
interactions occurring at high redshift, when the gas content 
is very high. In this sense, the merger remnant represents 
objects at z ~ 2 — 3 and the galaxies still have time to form 
more stars until z — 0. The subsequent infall of extragalactic 
gas ma y also produce more stars at z — and help to reform 
a disc ([Brooks et ai1l2009l ). 

3.1 Structure and Density Profiles 

At the end of each simulation, we calculated the V-band 
luminosity of the merger remnan t , usin g the stellar synthe- 
sis model of iKodama &; Arimotol (|l997h . From this, we pro- 
duced mock V-band luminosity maps (modulo the lack of 
dust extinction). Figs. [2Hl]show the luminosity maps for 3 
remnants. On each figure, the left panels show the stars born 
after the merger (young stars), the middle panels show the 
stars born before or during the merger (old stars) and the 
right panels show all stars. 

We immediately see that old and young stars have very 
different distributions. In simulation M12 (Fig. [2}, young 
and old stars form discs that have comparable radii, but the 
edge-on views (bottom panels) clearly shows that the young 
disc is thin, while the old disc is thick. The presence of two 
distinct discs, a thin one and a thick one, is in remarkabl e 
agreement with observations ([Yoachim &; Dalcantonll2006h , 
and this motivated us to publish the results of that partic- 
ular simulation in an earlier paper (B07). We found that all 
the simulations except for simulations M12z and Mil do re- 
sult in the formation of a thin disc made of young stars, and 
a thick disc made of old stars that were already present be- 
fore the merger, or formed during the mergers. These stars 
end up either in the thick disc or in the halo. 

Simulations M12z and Mil show more elliptical like 
remnants. Simulation M12z (Fig. [3} produced a remnant 
that has a very complex structure. Even though we find a 
small disc made of young stars, the overall structure resem- 
bles more an elliptical galaxy than a disc galaxy. Simulation 
Mil also shows a small young disc. Compared with the other 
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Figure 2. V-band image of M12 seen face-on (top panels) and edge-on (bottom panels), for young stars (left panels), old stars (middle 
panels), and all stars (right panels). X, Y, Z represent cartesian coordinates, with Z along the axis of rotation. 
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Figure 3. V-band image of M12z seen face-on (top panels) and edge-on (bottom panels), for young stars (left panels), old stars (middle 
panels), and all stars (right panels). The gas shows a complex structure that lasts for several hundred million years. 



simulations, the initial galaxies Gall and Gal2 are less mas- 
sive in simulation Mil, but the intensity of the starburst 
taking place during the merger is comparable. As a result, 
there is less gas available after the merger to form young 
stars and build up an extended thin disc. 

In most simulations, the merger resulted in the forma- 
tion of a ring made of young stars. This implies that pro- 
grade gas-rich major mergers have a tendency to form rings. 
These rings seem to be signatures of a gas substructure since 



most of the rings are coplanar with the thin disc. The large 
impact parameter use here is quite contrary to the normal 
scenarios forming a ring remnant ([Hernquist fc Weilll993l : 
iMapelli et al] |20Q8). That can be a signature of a gas-rich 
merger, and more work is needed to certify this assertion. 
The simulation M12z formed a very complex structure, with 
a non-coplanar ring connected to a very small thin disc by 
spiral arms, embedded in an ellipsoid distribution of old 
stars (Fig. [3}. Also, the face-on views of the mergers for 
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Figure 4. V-band image of Ml 290 seen face-on (top panels) and edge-on (bottom panels), for young stars (left panels), old stars (middle 
panels), and all stars (right panels). For this simulation, the stellar ring and the disc are not coplanar. 



to 

CM 



\ o - 



J* 

v — ' 
Q 



— i — i — i — | — i — i — i — i — | — i — i — i — r~_ 

M1290: 

KB ; 
r \^ i 

\ ^1 


71 — i — i — i — | — i — i — i — i — | — i — i — i — r~_ 

I rM12 j 

7 O*-.^ 

■ ° * * ; 
" i i i i I i i i i I i i i i " 


7T — i — i — i — | — i — i — i — i — | — i — i — i — r~_ 

I M12orbj 

: □ : 

- 

■ -.o-s : 

: e v.^j 

• * o. : 


— i — i — i — | — i — i — i — i — | — i — i — i — r~_ 

p M12 j 

" s - 

s. e : 

7 N^, Q □ ~ 

: □ : 

r ^ i 

: □ \ * 

" i i i i I i i i i I i T i I s 


; i i i i | i i i i | i i i i ; 

^. M12z : 
:o - 

7 OqO. " 
I ^ 

: * * . o 1 

Lj i i i 1 i i i i 1 i i @ .i 


: 
■_ M13 j 

: *V i 

1 * □ □ j 
: 

Lj i i i 1 i i i i 1 i i 


j i i i | i i i i | i i i r; 

b M110 : 

7 " 

ic2 : 

= ^%o ; 

r 

Lj i i i 1 i i i i 1 i i i i_I 


d ' ' : 
£ M11 j 

: □ : 

K i 

F ; 
- - 

: •<^. : 
: i i n 

1 1 1 1 1 1 1 1 1 1 1 1 ' " 



10 



10 5 
Radius (Kpc) 



10 



10 



Figure 5. Radial luminosity profile in the V-band, for old stars (squares) and young stars (circles). The straight lines show exponential 
fits of the old stars (dashed) and young stars (dotted). The dotted curves for simulation M12z and Mil are fits to a de Vaucouleurs 
profile for the young stars. 



simulations M1290 and rM12 (top panels of Figs. S] for an 
example ) clearly reveal the presence of a central bar made 
of both old and young stars. No bars are seen in the other 
simulations. 

We can use the surface density brightness to trace the 
radial luminosity profiles of the discs. The radial luminosity 
profiles are shown in Fig. [5] for all simulations. We calcu- 
lated these profiles separately for young stars and old stars. 
The dotted and dashed lines are fits to the profiles of young 
and old stars, respectively. In cases where a ring was clearly 
visible, we excluded it from the fit. All profiles are well- fitted 
by an exponential profile, except for simulations M12z and 



Mil, where the young stars are better- fitted by a de Vau- 
couleurs profile. We can use the exponential profiles to cal- 
culate the scale- lengths of the discs. In Table [3] we list the 
scale-length h of the old population, the scale-length h y of 
the young population, and their ratio h /h y , except for the 
simulations M12z and Mil, for which an exponential profile 
does not fit the distribution of young stars. The scale-length 
is usually larger for the old population than the young one, 
with ratios varying from 0.98 to 1.6 0. This is in agreement 
with t he recent observational work to lYoachim fc Dalcantonl 
(2006) if we make a direct correspondence between old stars 
and thick- disc stars. 



8 S. Richard et al. 

Table 3. Scale-length ratio for the young and old populations 



Run 


h [kpc] 


h y [kpc] 


h /K 


M12 


5.56 


5.27 


1.06 


M1290 


6.22 


6.05 


1.03 


M12orb 


8.97 


5.59 


1.60 


M12z 


6.55 






rM12 


7.54 


5.21 


1.45 


Mil 


5.94 






M13 


5.37 


4.47 


1.20 


MHO 


7.71 


7.87 


0.98 



Table 4. Angle between the plane of the disc for by the young 
and old populations 



Run 


6 oy [degrees] 


o [degree] 


&oy/ 0o 


M12 


3.28 


17.62 


0.19 


M1290 


0.70 


11.19 


0.06 


M12orb 


3.26 


20.64 


0.16 


M12z 


6.00 


41.17 


0.78 


rM12 


2.21 


1.27 


1.74 


Mil 


2.31 


26.68 


0.08 


M13 


5.63 


17.82 


0.32 


MHO 


0.25 


13.80 


0.02 



By definition, by the time the merger is completed, old 
stars have all formed, while the matter destined to form 
young stars is still in the form of gas. Hydrodynamical pro- 
cesses, such as oblique shock waves, can potentially exert a 
torque on the gas, modifying its angular momentum. Such 
processes would not affect the old population directly, only 
indirectly through the gravitational interaction between the 
old stars and the remaining gas. This might result in a mis- 
alignment between the young and old discs. To investigate 
this issue, we performed a least-square fit of a plane to the 
final distribution of young stars, and to the final distribution 
of old stars. The second column of Table 2] shows the angle 
O y in degrees between the two planes, that is, the angle be- 
tween the young and old discs. The angles are 6° or less for 
all the runs, indicating that the discs tend to be coplanar. 
To estimate whether these angles are significant, we fitted 
the distribution of stars in the old disc using an oblate ellip- 
soid. We then calculated a relevance angle o , defined as the 
arc tangent of the short-to- long axis ratio of the ellipsoid. If 
we were to fit the old disc inside a flat square box, o would 
be the angle between the edge and the corner, as seen from 
the center, so it is a measure of how thick the old disc is. If 
oy < o , then the young disc is entirely embedded inside the 
old disc, in spite of the angle between them. But if oy > o , 
then the young disc "sticks out" of the old disc. The third 
column of Table U shows the relevance angles. Interestingly, 
the simulation rM12, the only one for which one galaxy is 
spinning in the retrograde direction, produces an old disc 
that is actually over an order of magnitude thinner than the 
ones produced by the other simulations It is well-known that 
a retrograde disc will suffer very little ti dal disruption prior 
to m erger compared to a prograde disc ([Toomre &; Toomrd 
1972). This could explain why the old disc in simulation 
rM12 ends-up being quite thin, though a more detailed in- 
vestigation is needed to confirm this hypothesis. The last 
column of Table |4] shows the ratio oy /0 o . We find that all 
cases except rM12 have oy /0 o < 1. Hence, young discs tend 
to be embedded inside old discs. 

3.2 Kinematics 

We have analyzed the kinematic properties of all 8 remnants. 
We computed the rotation curves of the old and young popu- 
lations, for all the remnants. The top right panels of Figs. [6]- 
[9] show some examples of these curves. For simulations M12, 
M12orb, M1290, rM12, M13, and MHO, the old populations 
have a lower mean rotational velocity than their younger 
counterparts. The difference is about 200 — 300kms _1 for 



M12, M12orb, and M1290 (mass ratio 2:1), 150 kms -1 for 
M13 (mass ratio 3:1) and 50kms _1 for MHO (mass ratio 
10:1). Hence, it decreases with increasing mass ratio of the 
progenitors. Comparing M12 with rM12, we also find that re- 
versing the rotation of Gal2 (so that it is retrograde relative 
to the orbit of the galaxies) also reduces the differences be- 
tween the two populations (about 300 kins" 1 for M12 versus 
150 kms -1 for rM12). The rotation curve of M12z and Mil 
are drastically different. The young population of the outer 
ring is counter-rotating with respect to the young popula- 
tion of the outer disc. We also found some counter-rotating 
stars in other simulations, but (i) they belong to the old 
population, and (ii) these stars were in the minority, so the 
overall rotation curve showed no counter-rotation. 

To estimate the relative importance of rotation and ve- 
locity dispersion in providing support against gravity, we 
calculated for each population the average rotational ve- 
locity V and its dispersion a v in radial bins. The middle 
left panels of Figs. [BHS] the quantity ((V/ay) 2 ) 1 ^ 2 versus 
radius, for some of the remnants. In all simulations, we 
find ((V/a v ) 2 ) 1/2 > 1 for the young population, indicat- 
ing that this population forms a massive, rotationally sup- 
ported disc. The only exception is simulation M12z, where 
({V / cry) 2 ) 1 / 2 < 1 between radii 3 and 6 kpc. As the top 
panel of Fig. [71 shows, this region corresponds to the dis- 
continuity between the inner disc and the counter-rotating 
outer ring. Simulations M12, M12orb, M12z, and Mil all 
have {{V/ cry) 2 ) 1 / 2 < 1, indicating the presence of a "hot 
disc" supported by internal motions rather than rotation. 
Simulation M13 has ((V/a v ) 2 ) 1/2 < 1 at radii r < 3 kpc 
and ((V/a v ) 2 ) 1/2 > 1 at larger radii. Simulations M1290, 
rM12, and MHO all have {(V/a v ) 2 ) 1/2 > 1, indicating that 
the old population is also rotationally supported, but in all 
three cases the rotational support is significantly smaller for 
the old population, typically by a factor of order 5. In the 
case of simulation rM12, this result is consistent with the 
fact that retrograde orbits are less destructive than prograde 
orbits. 

The middle right panels of Figs. [BHS] show the veloc- 
ity dispersion versus formation epoch, with the dash verti- 
cal lines indicating the starburst. In 5 of the 8 simulations, 
the velocity dispersion is larger for stars formed prior to 
the collision, and smaller for stars born after the collision. 
This is consistent with the statement that old stars have 
a smaller rotational support than young stars. The excep- 
tions are the simulations with high mass ratios, M13 and 
MHO, for which a v peaks during the starburst, and M12z, 
for which a v peaks after the starburst (Fig. 0). 
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Figure 6. Star formation rate and kinematics for simulation M12. 
On each panel, the red curves represent old stars, the blue curves 
and dots represent young stars, and the black curves represents all 
stars. The vertical dashed lines indicate the beginning and the end 
of the merger. Top left: star formation rate vs. formation epoch. 
t = corresponds the the beginning of the simulation. Top right: 
rotation velocity vs. radius for stars located near the plane of the 
disc (\Z\ ^ lkpc). Middle left: Rotational support versus radius. 
The dashed line separates rotationally-supported stars (above) 
from stars supported by velocity dispersion (below). Middle right: 
velocity dispersion vs formation epoch. Bottom left: histogram of 
stellar mass vs. rotation velocity, withe negative values indicating 
counter-rotating stars. Bottom right: Toomre diagram. The black, 
green, and blue dots indicate stars formed before, during, and 
after the merger. 



The bottom left panels of Figs. [BHS] show histograms 
of the stellar mass versus rotation velocity. For simulations 
M12, M12orb, M1290, rM12, and M13, the young population 
is concentrated in a narrow region of the histogram, while 
the old population is spread in velocity, from the positive to 
the negative side (though it is concentrated mostly on the 
positive side for rM12). The young component is counter- 
rotating for simulation Mil. For simulations M12z, parts 
of the young component are counter-rotating, and the to- 
tal distribution is centered near zero. Simulation Ml 10 (the 
minor merger) differs from all others in that the old popu- 
lation is concentrated at positive velocities. Notice that this 
is the only case for which the old population is rotationally 
supported. 

One way to discriminate between kinematically differ- 
ent star populations i s the use of the Toomre Diagram 
(jSandage fe Fo uts 1987). By comparing V (circular velocity) 
with U-\-W (radial and perpendicular velocity respectively), 
we can have a clearer picture of the kinematic distribution 
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Figure 7. Same as Fig. [6] for simulation M12z. 
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Figure 8. Same as Fig. [6}{7] for simulation rM12. 
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Figure 10. Circularity of the orbits for all simulations, ej rep- 
resent the angular momentum of a star particle normalized to a 
particle in a circular orbit of the same radius. On each panel, the 
red curves represent old stars, the blue curves represent young 
stars, and the black curves represent all stars. 



Figure 9. Same as Fig. [6j|HI f° r simulation M13. 



of the young and old population. A thin disc will be mostly 
formed by stars with low U + W and a large V compo- 
nent since all stars will be in nearly coplanar orbits. The 
bottom right panels of Figs. [6H2]show the Toomre diagram 
for each simulations. For M12, the young stars are concen- 
trated in a small region of the diagram, centered around 
V = 400 km s" 1 and (U 2 + W 2 ) 1/2 = Okms" 1 , which is 
characteristic of a thin disc since these stars are mainly 
on fast rotating planar orbits. The old stars are distributed 
throughout the diagram, indicating that these stars follow 
orbits with significant radial and orthogonal velocities, char- 
acteristic of systems supported by velocity dispersion. These 
stars are located mostly in the thick disc and halo. 

The Toomre diagrams for the other simulations show 
similarities, and also interesting differences, with the dia- 
gram for simulation M12. Some simulations present a dis- 
continuity between the positive and negative region of the di- 
agram. For example, simulation rM12 produces two different 
population of old stars, one forming a thick disc and another 
one forming a counter- rotating thick disc/halo. Simulation 
M12z and Mil have a different Toomre diagram than the 
other simulations. For M12z (Fig. 0), there is no clear disc 
formed in the remnant, and the young stars are dispersed 
throughout the diagram. As we saw in the previous section, 
young stars in this simulation are mostly located in a mas- 
sive ring around the galaxy. For simulation Mil, the young 
stars in the outer ring are counter-rotating while all stars 
in the inner disc are co-rotating. For simulation Ml 10, the 
young population has a significant (U 2 + W 2 ) 1 ^ 2 component 
(rsj 150 kms -1 ), almost as large as the V component. 

We have also performed an analysis of the circularity 
of the stellar orbits, where ej = j z /j c irc, and j z is the z- 



component of the specific angular momentum of each star, 
and j c i rc is the angular mom entum expected for a circular or - 
bit, following the method of lSanchez-Blazquez et all (12 009), 
which is an approximation of a method used bv lAbadi et al.l 
(2003). These results are showed in Fig. [TUJ Simulations 
M12, M13, MHO, M1290, rM12, and M12orb all present 
a young star population with circular orbits (ej ~ 1.0) and 
an old population with a high dispersion of circularity. One 
exception is MHO, which presents near circular orbits for 
both the young and old population, which is expected in a 
minor merger. Simulations Mil and M12z show non-circular 
orbits, in good agreement with the other kinematical re- 
sults. The small thin disc of Mil is also well represented at 
e, ~ -1.0. 

The main conclusion we can draw from the analysis and 
comparison of the kinematical properties of the remnants 
is that these kinematical properties, as the structural ones, 
are strongly depen dent of the initial c ondit ions, and support 
previous results bv lRobertson et al.l (|2006h . 

3.3 Chemical Abundances 

As mentioned in the Section I, one of the most important 
results not yet achieved concerning gas-rich major mergers 
is their chemical signature. Here we present some of our 
first results. The chemical properties of all remnants are 
presented in Figs. [TTHT81 The top panels show the radial 
profiles of the a-element abundances [a/Fe] (calculated by 
averaging the abundances of O, Mg, and Si) and metallicity 
[Fe/H] , calculated by averaging over radial bins. The old 
population has a higher [a/Fe] ratio and lower metallicity 
[Fe/H] than the young population, at all radii. 

The gradients of [a/Fe] in the old population are all very 
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flat. The distributions of [a/Fe] in the young population are 
much more complex. In particular, simulations Ml 2, M12z, 
M1290, and M13 show several local minima and maxima 
at various radii, with variations as large as 0.5 dex. The 
gradients of [Fe/H] are much smoother. In all simulations 
except Mil, we find a gradient in [Fe/H], ranging from —0.01 
to —0.1 dex per kpc. 

The third rows of panels in Figs. [TTHT8l show the [a/Fe] 
and [Fe/H] ratio versus formation epoch. The [a/Fe] ratio re- 
mains fairly constant over time during the merger, and starts 
decreasing after the merger. The [Fe/H] ratio shows a com- 
pletely opposite behavior: it increases strongly before and 
during the merger, and keeps increasing, but more slowly, 
after the merger. This would suggest that the galaxy grows 
outside-in, with the older, metal-poor stars being located at 
the outer regions of the disc. To investigate this question, we 
plot in Fig. [19] the formation epoch of old and young stars 
versus their final location in the disc, for all simulations. 
The dashed lines indicate the beginning and the end of the 
merger. For old stars, the formation epoch either decreases 
with increasing radius (as in Ml 2) or remains constant (as in 
M1290). For instance, in the case of M12, stars that formed 
before the merger are dominant at radius R>3 kpc, while 
stars that formed during the merger and the associated star- 
burst are more centrally concentrated. The young stars show 
the opposite trend: the formation epoch increases with ra- 
dius up to R = 7.5 kpc in most simulations. Star-formation 
after the merger therefore proceeded inside-out. It is the ef- 
ficiency of chemical enrichment, and not the epoch of star 
formation, that explains the gradient in [Fe/H]. Chemical 
enrichment is more efficient in the center of the remnant, 
where the stellar density is higher. 

The second row of panels in Figs. [TTHTHl show [a/Fe] 
and [Fe/H] versus height above the plane. The different ex- 
tents of the two curves on each panel reflects the different 
thicknesses of the two discs. Notice that the thin disc is very 
thin in simulation Mil, less than 800 pc. There is essentially 
no vertical gradient in [a/Fe] for the old population. For the 
young population, we find a gradient in [a/Fe] for M12 (0.1 
dex per kpc), and possibly another one for M12orb, while 
there are no evidence for gradients in the other cases. The 
[Fe/H] ratios are either constant, as in simulations M12z and 
Ml 10, or drops by a factor of order 2 or 3, as in simulations 
M12 and M1290, for both the old and the young populations. 
Since the stellar density, and the resulting chemical enrich- 
ment, is a much stronger function of radius than height, we 
naturally expect the vertical gradien t of the metallicit y to 
be much weaker than the radial one (jBrook et alJl2005h . 

The most important plots on Figs. Illf[18l are the [a/Fe] 
versus [Fe/H] plots (bottom- left panels). Old stars have 
larger a-element abundance than young stars, up to a rel- 
atively high metallicity ([Fe/H] ~ —0.5). We can explain 
this result by considering the relations between [a/Fe], 
[Fe/H], and the formation epoch of stars (third- row panels 
in Fig. I11H18[) . The metallicity [Fe/H] increases with time 
during the starburst, and levels-off after the starburst is 
completed. The ratio [a/Fe] then decreases with time as the 
metallicity increases, until it reaches the values found in the 
thin disc. This decrease of [a/Fe] is slowed near the begin- 
ning of the collision, as the starburst leads to a large number 
of Type II SNe, that enrich the gas in a-elements. After the 
starburst, Type la SNe become effective, and enrich the gas 
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Figure 11. Chemical abundances for simulation M12. On the 
top four panels and bottom left panel, solid curves represent 
old stars, and dot-dashed curves represent young stars. The first 
three rows show the a-elements abundance (left) and metallicity 
(right) versus radius (first row), height (second row), and forma- 
tion epoch (third row). Bottom left panel: ai-element abundance 
versus metallicity. Bottom right panel: Stellar mass versus metal- 
licity, for stars born before or after the merger (medium line), 
after the merger (thin line), and for all stars (thick line). 



in iron. This explains why the ratio [a/Fe] gradually de- 
creases after the collision. 

We complete this study by calculating the stellar mass 
in metallicity bins, for both populations. The results are 
shown in the bottom right panels of Figs. 11114181 Old and 
young stars have different metallicity distributions. The 
old population peaks in the range [Fe/H] = [—1.0, —0.6], 
while the young population peaks in the range [Fe/H] = 
[—0.4,0.4]. The implies that old stars are formed, on av- 
erage, by gas that has not yet been enriched by a large 
number of nearby Type la SNe. We can gain more insight 
into this process by plotting Toomre diagrams for different 
metallicity bins. The results for simulation M12 are shown 
in Fig. 1201 We recall that the concentration of stars in the re- 
gion V ~ 400km/s, (U 2 + W 2 ) 1/2 ~ Okm/s represents the 
thin disc, and the reminder of the diagram represents the 
thick disc and the halo. Low- metallicity stars ([Fe/H] < —3) 
are located in the halo, high-metallicity stars ([Fe/H] > 0) 
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Figure 12. Same as Fig. [TT] for simulation M12orb 
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Figure 13. Same as Fig. [TTlfHl for simulation M12z 



are located in the thin disc, and intermediate-metallicity 
stars are found in both locations. 

The main result we draw from this analysis of the chem- 
ical abundance is that the merger leaves a clear signature: 
the old stars, located in the thick disc and the halo, have a 
ratio [a/Fe] that remains constant with increasing [Fe/H] up 
to high metallicities ([Fe/H] = —0.5), while the young stars, 
located in the thin disc, have a lower ratio [a/Fe] which 
decreases with increasing metallicities. This result is very 
robust: we found it in each of the 8 simulations included in 
our study, including the minor merger (simulation Ml 10). 

Early gas-rich mergers in our simulations occur, by con- 
struction, prior to the timescale for significant Type la SNe 
pollution. Ultimately, this is responsible for the robustness 
of our results. 



4 DISCUSSION 

The results presented in this paper support the gas-rich 
merger scenario as a way to explain the formation of large 
disc galaxies in the Universe . Our results agree ver y well 
with the kinematical results of lRobertson et al.l ([20061 ) . That 
said, the gas-rich process may be necessary to form a disc, 
but it is not a sufficient condition. We can see this in our two 
simulations which formed bulge-dominated galaxies, namely 
Mil and M12z. Explaining this difference is not simple be- 



cause of the chaotic nature of the merger process. There is 
likely a complex relationship between the initial conditions 
and the behavior of the gas during the merger. The dura- 
tion of the merger and the violence of the collision seems 
to favor a higher dynamical warming, and this favors the 
formation of a small rotationally supported disc, embedded 
within a large spheroidal component. In the case of simu- 
lation Mil, the two progenitors have a lower mass than in 
the others simulations, but the starburst has the same in- 
tensity. As a result, there is not enough gas left after the 
merger to form a substantial thi n disc. This agrees with the 
results of lRobertson et all (|2006l ) who found that low-mass 
progenitors are more likely to form an important spheroidal 
component than high-mass ones. 

When a disc- like remnant is formed, it invariably con- 
sists of two disc components: a thick disc made of stars 
formed before or during the merger, and a thin disc made 
of stars born after the merger. The thick disc has a lower 
rotation ve locity than the thin disc, in a greement with the 
findings of lYoachim fe Dalcantonl ([2006). The scale-length 
ratio between the two components is also in agreement with 
this work, which suggests that gas-rich disc mergers may be 
a common way to form disc galaxies. However, we have to be 
careful since subsequent gas infall from satellites or the in- 
tergalactic medium is not included in our simulations. Such 
gas infall could increase t he scale- length of th e thin disc and 
affect our results (see, e.g. jBrooks et al . 2009). Stars belong- 
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Figure 14. Same as Fig. [TTlfT3l for simulation M1290. 
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Figure 15. Same as Fig. [TTlfT4l for simulation rM12. 



ing to the thick and thin discs differ by their abundances in 
various metals. Old stars have an excess in a-elements rela- 
tive to young stars at metallicities [Fe/H] below —0.5. This 
result is in good agreement with observational studies on 
the chemical abundances of the two disc components of the 
Milky Way, if the old and young stars in our simulations 
are analogous to the thick- and thin-disc stars, respectively. 
iReddv. Lambert, fe Allende Prietol (|2006h published a spec- 
troscopic survey of 176 stars with high probability to be 
part of the galactic thick disc, selected with the Hipparcos 
catalog. These authors analyze the relationship between the 
abundances of 22 chemical elements and the [Fe/H] ratio. 
It is clear form their Fig. 12 that the [a/Fe] ratio is higher 
for the thick disc than the thin disc by almost 0.15 until 
the metallicity reaches [Fe/H] = —0.3. At larger metallici- 
ties, the [a/Fe] ratios of the two discs becomes similar. The 
authors explain this fact by invoking an extended period of 
enrichment in iron by Type la SNe, an explanation that is 
supported by our work. Fig. 16 in the same article shows 
that chemical elements other than the a-elements do not 
follow the same pattern, which is further evidence of the 
important role of Type II SNe in the formation of the Milky 
Way, and support the scenario that a major collision took 

place in the history of the Milky Way. 

Another similar study was done bv lBensbv et al.l f|2Q05h . 
That study presents abundances for 102 dwarf F and G 
stars, and uses their kinematics to determine if they belong 



to the thin or thick disc. Their Fig. 8 shows the abundances 
of oxygen, magnesium, and silicon as functions of metallic- 
ity, which indicates that that stars in the thick and thin disc 
have different abundances for a metallicity [Fe/H] < —0.5. 
Their Fig. 10 is particularly interesting and shows that the 
abundance of oxygen decreases with increasing metallicity 
in the interval —1.0 < [Fe/H] < 0.0. This is all in agreement 
with the results presented here. These authors suggest the 
existence of several observational constraints related to the 
formation and chemical evolution of the thick disc: 

• Thick-disc stars and thin-disc stars have different chem- 
ical abundances. 

• For a [Fe/H] smaller than a certain value, thick disc 
stars have a larger abundance of a-elements than thin-disc 
stars. 

• The ratio [a/Fe] decreases with increasing [Fe/H], 
which shows the important contribution of Type la SNe. 

• Stars in the thick disc are on average older than stars 
in the thin disc. 

Our results mostly agree with these constraints, which 
supports the gas-rich, major- collision scenario for the forma- 
tion of the thick disc of the Milky Way, a scenario in which 
a collision leads to an intense starburst accompanied by a 
very large number of Type II SNe. These SNe enrich the gas 
supply in a-elements, allowing the formation of stars with a 
high [a/Fe] ratio. The starburst also leads to a progressive 
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Figure 16. Same as Fig. fTTTfTBI for simulation Mil. We could 
not calculate vertical profiles for the young stars, because the thin 
disc is too thin. Instead we indicate the mean values of [a/Fe] and 
[Fe/H] with a single symbol (x) in the second-row panels. 



enrichment in iron caused by Type la SNe. Stars formed af- 
ter the collision will therefore have a smaller [a/Fe] for the 
same metallicity. These stars will form a thin disc with a 
smaller velocity dispersion than the stars formed during the 
merger. Another interesting result is the drop in [Fe/H] with 
increasing radius. 

Finally, we must point out that the initial conditions 
used in our study are greatly simplified. First, the initial 
conditions do not include a spheroidal component made 
of stars. We assume that the merger progenitors are pure 
disc for simplicity. We might expect that these stars, if in- 
cluded, would end up in the halo of the final galaxy, which 
might explain that most simulations produced a low-mass 
halo. We can argue that our results are valid in the limit 
where the initial spheroidal components are too small to 
affect the dynamics of the collision. We have not included 
a bulge component as well. Although morphologies of high 
redshift galaxies are unknown, this assumption may be too 
simplistic and our results may be biased by this assump- 
tion. However, a recent work showed the presence of a bulge 
can produce a disc remnant with a gas ratio as low as 12% 
(|Barnesl l2002). These facts are an evidence that a merger 
between bulgeless disc-disc galaxies is not a sufficient con- 
dition to form a dick remnant. In fact, bulgeless mergers 
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Figure 17. Same as Fig. [TTlfT6l for simulation M13 



generall y produce elliptical remnants in di ssipationless simu- 
lations ([Gonzalez- Garcia &; Balc ells 2005). So, the presence 
of a massive gaseous component is probably the major con- 
dition to reform a disc after a major merger. We also have 
to add that the softening length used in this study is slightly 
higher than what are used in the recent state-of-the-art sim- 
ulations. However, the morphological aspects of our study 
show good agreement with previous such higher spatial res- 
olution simulations, showing that the resolution is sufficient 
for the analysis presented in this paper. 

The hierarchical model of structure formation normally 
includes a significant amount of collision and accretion of 
low-mass galaxies. It is presumptuous to assume that the 
formation of the thick disc can be explained entirely by one 
single collision. It is possible that several of these collisions 
take plac e in the initial phases of the form ation of spiral 
galaxies (Bro ok et al1l2005l : IConselicd l2006h . Nevertheless, 
our study can explain several observations with the impor- 
tant starburst that takes place during a major, gas-rich col- 
lision. 



5 SUMMARY AND CONCLUSION 

Using GCD+, we have performed 8 simulations of major 
mergers between gas-rich spiral galaxies. We have analyzed 
the kinematic, structural, and chemical properties of stars 
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formed before and during the collision (the old population) 
and stars formed after the collision (the young population). 
We used the star formation rate to define these two pop- 
ulations. A fraction of the old stars end up in the halo of 
the merger remnant, while the remaining stars form a thick 
disc which is partly supported by velocity dispersion, partly 
supported by rotation, and sometimes includes a significant 
counter-rotating component. The young stars form a thin 
disc that is supported by rotation, and, in many cases, a 
ring that might or might not be coplanar with the thin disc. 
The disc themselves tend to be coplanar, the angle between 
them being ~ 6° or less. With rare exceptions, both discs are 
well-fitted by an exponential profile, and the scale-length of 
the thick disc exceeds the one of the thin disc a few percents 
up to a factor of 1.60. 

The starburst occurring during the collision rapidly en- 
riches the gas in various metals. Explosions of Type II SNe 
follow rapidly the start of the collision, owing to the short 
lifetime of their progenitors, and enrich the intergalactic 
medium in a-elements. Enrichment by Type la SNe is spread 
over a large period of time, enabling the enrichment in iron 
of both stellar populations, which results in an old popula- 
tion having a [a/Fe] ratio higher than the young population, 
even at relatively large metallicities ([Fe/H] = —0.5). This 
could explain the high [a/Fe] ratio observed in stars in the 
thick disc and the halo of the Milky Way. This result do 



not depend strongly upon the initial conditions, since it was 
found in all simulations. 

Our main conclusion is that the morphological, kine- 
matical and chemical properties of the thick and thin disc 
can be reproduced in an scenario where thick disc formed in 
a gas rich merger between two disc galaxies. Furthermore, 
while the structural and kinematical properties of the merger 
remnants are strongly dependent on the initial conditions - 
features such as the ratio of the disc scale- lengths, the extent 
of the discs, the presence of rings or of a counter-rotating 
component, vary from simulation to simulation- the chem- 
ical abundances show a remarkable consistency among the 
various simulations. The key result of this study is that the 
ratio [a/Fe] remains constant with increasing [Fe/H] for old 
stars, up to [Fe/H] = —0.5 while it decreases with increas- 
ing metallicity for young stars. This is true in each of the 
8 simulations we considered. The observed chemical signa- 
tures in our merger remnants require that mergers happen 
before the onset of the majority of Type la SNe, regardless 
of the orbits or mass ratios of the progenitors. Furthermore, 
the mergers need to be gas rich to lead to a remnant with a 
rotationally supported disc. These characteristics are consis- 
tent with the mergers we expect to find at high redshift and, 
therefore, our results support this scenario for the formation 
of the different disc components in the Milky Way. 
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